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Abstract 

We  consider  Toda  flows  induced  on  the  set  of 
orthogonal  upper  Hessenberg  matrices.  The 
explicit  formulas  for  the  evolution  of  Schur  pa¬ 
rameters  are  given. 

1  Introduction 

Any  symmetric  nonnegative  definite  Toeplitz 
matrix  rn+i  of  order  n  +  1  can  be  modeled 
as  an  autocorrelation  matrix  of  a  stationary 
signal  [12] 

p 

xm  =  Yla‘  «>s (rnu>i  +  &i)  +  ym, 

/= l 

where  0/  are  arbitrary  phase  shifts  and  ym  is  a 
zero  mean  white  noise  process  whose  variance 
equals  the  smallest  eigenvalue  Amin  of  Tn+j. 
Assume  that  the  eigenvalue  A,„in  is  simple,  and 
let  (t/o,  •  •  • ,  fjn)  be  a  corresponding  eigenvector. 
Then  [12]  the  polynomial 

V>n(A)  =  »*)  +  ...+  »?„An 

has  n  distinct  roots  Aj,*  •  -  ,An  on  the  unit  cir¬ 
cle,  and  the  frequencies  of  xm  are  given  by 
{exp(±iu»/}{’_1  =  {Aj}”=1,  where  i  denotes  the 
imaginary  unit.  One  can  construct  [2]  an  or¬ 
thogonal  Hessenberg  matrix  0  with  character¬ 
istic  polynomial  proportional  to  V’n-  Moreover, 


the  amplitudes  a/  can  be  recovered  from  the 
first  components  of  the  normalized  eigenvec¬ 
tors  of  O.  One  can  then  use  any  of  several  al¬ 
gorithms  designed  for  unitary  and  orthogonal 
Hessenberg  eigenproblems  [9,  6,  1,  10,  4] 
to  calculate  the  frequencies  and  amplitudes. 
In  the  present  paper  we  investigate  another 
aspect  of  orthogonal  Hessenberg  matrices. 
Namely,  we  consider  Toda  flows  on  these  ma¬ 
trices  (referred  to  as  Schur  flows  in  [3])  and 
obtain  explicit  formulas  for  the  evolution  of 
the  so-called  Schur  parameters  under  the  Toda 
flow.  Since  Schur  parameters  determine  or¬ 
thogonal  Hessenberg  matrices  uniquely,  we  ac¬ 
tually  obtain  an  explicit  description  of  the  evo¬ 
lution  of  a  given  orthogonal  Hessenberg  matrix 
under  the  Toda  flow. 

2  Inverse  problem  for  orthog¬ 
onal  Hessenberg  matrices 

Let  M  be  the  set  of  positive  Borel  measures 
on  C  which  have  the  following  properties.  For 
any  n  €  M  the  support  of  /i,  which  we  denote 
by  Am,  consists  of  exactly  n  points  which  lie 
on  the  unit  circle  U  and  is  such  that:  i)  if 
A  6  Am,  then  the  complex  conjugate  A  is  also 
in  and  /t{A)  =  /x{A);  it)  jt{A)  >  0  for  any  A 
in  AM;  tit)  n(C)  =  1;  iv)  -10  AM.  We  further 
introduce  a  class  OH+  of  orthogonal  matrices 
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O  =||  Oij  ||  such  that  o,j  =  0  if  *  —  j  >  1  , 
o.+i,«  >  0  for  all  i  and  det  0  =  1.  Finally,  given 
a  vector  r  =  (r0,  •  •  • ,  rn_ i  )T  £  Rn  introduce  a 
corresponding  Toeplitz  matrix  T(r)  =||  t,j  ||, 
where  ttJ  =  T|,_j(. 

Theorem  2.1  Given  a  positive  definite  sym¬ 
metric  Toeplitz  matrix  T(r)  with  tq  —  1  there 
exist  exactly  one  measure  p  £  M  and  exactly 
one  0  £  0H+  such  that 

f  X*dp(X)  =  t,-  =<  e\,Ole\  >,  (2.1) 

Jc 

i  =  0,  •  •  • , n  —  1.  Here  e\, •  •  • , en  is  the  canoni¬ 
cal  basis  in  Rn  and  <,  >  is  the  standard  scalar 
product.  Conversely,  for  any  p  £  M  and  any 
O  £  OH+  the  matrices  T(r),  T(r')  are  Posi¬ 
tive  definite  Toeplitz  matrices.  Here 

t{  =  f  X'dp(X),rl  =<  euO'ei  >, 

J  c 

i  =  0,  •••,«  —  1. 

Remark  2.2  Theorem  2.1  is  more  or  less 
known  to  the  experts  (see  e.g.  [8],  [11]).  We 
nevertheless  give  an  independent  proof  to  clar¬ 
ify  relationships  between  introduced  objects. 

Remark  2.3  There  is  nothing  mysterious 
about  the  number  —1  which  we  have  excluded 
from  the  support  of  each  measure  in  M.  This 
simplifies  notations  a  little  bit. 

We  need  the  following  elementary  lemma. 

Lemma  2.4  Let  uj,* ••,vn-i  be  an  orthonor¬ 
mal  system  of  vectors  in  Rn.  There  exists  ex¬ 
actly  one  orthogonal  matrix  O  such  that  Oei  = 
«,-,*  =  1, •  ••  ,n  —  1  and  detO  =  1. 

We  can  now  outline  a  proof  of  Theorem  2.1. 
Proof:  Denote  by  Pn  the  vector  space  of  real 
polynomials  of  degree  less  or  equal  n  -  1.  Set 

<A*,A'>=  /  X'-idp(X).  (2.2) 
Jc 

We  prove  that  (2.2)  defines  a  positive  definite 
scalar  product  on  P„.  Observe  that 

J  X'X’dp  =  J  X {~jdp  =  J  XXjdp. 


Indeed,  /  X'X’dp  =  A*A^{A)  = 

A*-j/u{A},  since  A  =  A-1  for  A  £  U. 
Further,  since  /z{A}  =  p{ A)  we  have 

£  X'-ip{X}  =  £  A-MA}  =  /  A’ A  Up. 

Let  q  =  ao  +  . . .  +  an_iAn_1  £  Pn.  We  have 

<?,?>=  £  a«°i  /  A’AJd/i  =  [  |  q  |2  dp  >  0. 

*J= o  J  J 

Further,  /  |  q  |2  dp  =  0  if  and  only  if  q( X)  =  0 
for  any  A  €  AM.  Since  degq  <  n  =  card( AM), 
this  is  possible  only  if  q  =  0.  Consider  the 
polynomial  £(A)  =  n«eAM(A  -  t)  =  b0  +  . . .  + 
6„_iAn_1  +  A”.  Since  all  roots  of  £  lie  on 
the  unit  circle  we  clearly  have  A"f(l/A)  = 
&of(A),f>o  =  ±1-  Further,  all  coefficients  of 
£  are  real  because  AM  =  AM.  Consider  the 
linear  operator  O  :  Pn  —>  Pn  defined  as  fol¬ 
lows:  OA*  =  A,+1,i  =  0, •••,«  — 2,  OA”-1  = 
—bo~b1X  — . . .  —  hn_; A"-1 .  We  now  prove  that 
0  is  orthogonal  relative  to  the  scalar  product 
< ,  >  .  We  should  prove  that 

<  O X\X>  >=<  A',0_1AJ  > 

for  any  i,j  =  0,  •••,«—  1.  The  only  non¬ 
trivial  case  is  x  =  n  —  1,  j?  =  0.  We  have 
<  OA"-1, 1  >=  -b0  -  biTi  -  ...  -  fcn.jTn-!, 
where  ry  =  fc  X'dp(X).  Let  O-1 1  =  Co  +  . . .  + 
cn_,An-1.  Then  <  An-1,0-1l  >=  corn_i  + 
CitV»-2+-  •  -+cn_i.  Thus,  it  is  sufficient  to  prove 
that  a  =  — 6n_i_,,  i  =s  0,  •••,»—  1.  We  clearly 
have  1  ss  coOl  + . . .  +  cn_iOA"-1  =  coA  +  ...+ 
Cn-2An-1-t-c„_i(-6o-6iA-..  .-6n_iAn-1)  or 
1  =  -C„_i60,  Co  -  C„_i&i  =0,  C\  -  Cn_i&2  = 
0,  •  •  • ,  cn_2  -  cn_i6„_i  =  0.  This  yields  hi  = 
-co/ho,  h2  =  -ci/ho,  —  -,  hn_i  =  -cn_2/h0. 
We  now  use  the  relation  A"£(l/A)  =  ho£(A).  It 
follows  that  h„_,-  =  bobi,  i  =  0 ,•••,»».  Thus 
bn-i/b o  =  -Ci-i/bo,  i  =  1,*  •  *,n.  These 
are  exactly  the  required  conditions.  Thus 
we  have  constructed  an  orthogonal  operator 
0  such  that  fcX'dp  —<  1,0*1  >,  i  — 
0,  1.  Observe  that  the  characteristic 
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polynomial  of  O  coincides  with  f.  Thus  the 
spectrum  of  O  is  AM.  In  particular,  detO  = 
1  (here  we  use  the  assumption  that  —1  £ 
Am).  Let  po  =  1, •••,?„_!  be  an  orthonor¬ 
mal  basis  in  Pn  obtained  by  the  orthonor¬ 
malization  of  the  basis  1,  A,- •  ■ ,  A"-1.  It  is 
clear  that  the  matrix  O  of  the  operator  O 
is  upper  Hessenberg  in  this  basis.  Moreover, 
the  entries  di+ij  are  all  nonzero  (otherwise, 
span(po,  •••,p,_j)  =  span(l,-  •  -,A‘_1)  is  an 
invariant  subspace  of  O  which  is  not  true). 
Without  loss  of  generality  one  can  suppose 
that  d,+iit  >  0  for  all  i.  Otherwise  one  can 
take  diag(e j,  •  •  -,en)Odiag((i,  •••,€„). 

Suppose  we  are  given  a  positive  definite 
Toeplitz  matrix  T(r)  and  an  orthogonal  ma¬ 
trix  O  €  OH+  such  that  r,  =<  >,i  = 

0,  •  ■  -  ,n  —  1.  Then 

T(t)  =  VTV,  (2.3) 

where  V  is  the  upper  triangular  matrix 
[«i,Oei, •  ••  ,On_1ei]  with  positive  entries  on 
the  main  diagonal.  But  (2.3)  is  the  Cholesky 
decomposition  of  T(r).  Hence  it  is  uniquely 
defined  by  T(t).  In  other  words,  the  vectors 
Oei,’  •  -,On-1e1  are  uniquely  defined  by  T(r). 
Since  these  vectors  form  a  basis,  the  vectors 
Ot\ ,  •  •  • ,  Oen_i  are  uniquely  defined  by  our 
Toeplitz  matrix.  Thus  by  Lemma  2.4  the  ma¬ 
trix  0  is  uniquely  defined  by  T(t).  Given  a 
positive  definite  Toeplitz  matrix  T(r)  we  can 
endow  Pn  with  a  scalar  product  <,  >  and  the 
shift  operator  defined  on  span(l,  A,  ••  • ,  A”-2) 
as  we  did  before.  Then  using  Lemma  2.4  we 
can  extend  this  operator  to  the  orthogonal  op¬ 
erator  0,  defined  on  Pn  such  that  detO  =  1. 
Then  the  matrix  of  O  in  the  basis  obtained  by 
orthonormalization  of  the  basis  1,  A,-- -.A"-1 
belongs  to  OH+  and  77  =<  ei.O'ei  >,t  = 
0,  •••,«  -  1.  Consider  now  the  rational  func¬ 
tion 

/(*)=<  1»(*7  -  0)~l\  >  . 

As  is  easily  seen 


where  all  r,  >  0.  We  then  can  define  the  mea¬ 
sure  p  e  M  by  the  conditions  p{A,}  =  t\  and 
equal  to  zero  otherwise.  We  immediately  see 
that  equations  (2.1)  are  satisfied.  It  remains  to 
prove  that  the  measure  p  is  defined  uniquely 
by  conditions  (2.1).  Let  pk  G  M,k  =  1,2  be 
such  that 

[  X'dpi  =  f  \'dp2, 

Jc  Jc 

i  =  0,  •••,n  —  1.  Then  we  can  construct 
Ok,k  =  1,2  such  that  conditions  (2.1)  are 
satisfied.  But  then  Oi  =  O2.  In  particular, 
AMi  =  Am,  i.e.,  pi  =  p2  because  we  have 
for  p{A}  the  following  system  of  Vandermonde 
equations: 

A*/i{A}  =  r„ 

*  =  0,  ••*,»—  1.  ■ 

Let  T(t)  be  a  positive  definite  n  x  n  Toeplitz 
matrix  and  <,>  be  the  corresponding  scalar 
product  on  Pn.  Let 

p,( A)  =  A’  +  . . . ,  Si  >  0,  i  =  0,  •  •  • ,  n  -  1, 

be  the  basis  obtained  by  the  orthonormaliza¬ 
tion  procedure  from  the  basis  1,A, ••  •, An_1. 
Since  p, 

is  orthogonal  to  span{  1,  A,  •  •  • ,  A’-1),  we  have: 
Ap,(A)  is  orthogonal  to  span( A,  •  -  * ,  A*).  Fur¬ 
ther,  r  =  Ap,(A)/£,  -  pi+i/Si+i  6  Pi+i-  Let 

<Pi  G  Pi+\  be  such  that  <  q,(pi  >=  q(0)  for 
any  q  £  Pi+i.  Since  p,  is  orthogonal  to  Pi  ard 
both  r  and  <p,  are  orthogonal  to  APt,  we  obtain 

Ap,(A)/£,  =  pi+1(A)/£,+1  -»-  7,<p„  (2.4) 

for  some  real  7,,  *  =  0,  •  •  • ,  n  —  2.  An  easy  cal¬ 
culation  shows  that  <p,  =  0,A'p,(l/A).  Hence 

1  =  *</«*! +7  (2-5) 

In  other  words,  if  we  know  70, •••,7n-25  we 
can  find  6\ ,  •  •  • ,  6n-i  •  Then  using  (2.4),  one  can 
determine  p\,  •  •  •  ,pn-i  and  consequently  using 
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again  (2.4)  the  corresponding  upper  Hessen- 
berg  orthogonal  matrix  O.  We  have  by  (2.4) 
<  *?.(*)> P.( A)  >=  7Ap»(0),*  =  0,  -••,»-  2. 
Evaluating  (2.4)  at  0,  we  obtain  p,+i(0)  = 
-7«A?$*+ i*  *  =  0,  •  •  • ,  n  -  2.  Thus  o<+i,,+i  =< 
Ap,A,p,(A)  >=  «  =  l,---,n-2. 

Further,  oj, i  =  -7oPo(0)  =  ~7o-  Let  us  set 

a\  —  =:  ^t— l  /^*  ■»  v\  ~  7»-i^f-i» 

*  =  1,  •  •  • ,  n  —  1.  We  obviously  have 

a]  +  v}  =  1,  Oij  =  -I'.-ii /,, 


uniquely  defined  by  its  auxiliary  Schur  param¬ 
eters  <r,(t)  =  Oj+i,,-(<).  We  now  describe  ex¬ 
plicitly  how  these  parameters  evolve  under  the 
Toda  flow. 

Theorem  3.1 


i  =  1,  •••,»—  1,  Ao  =  1.  Here  A ;(<)  is 
the  i—th  principal  minor  of  the  matrix  T(t)  = 
exp((O(0)  +  O(0)T)t). 


*  =  n  -  1,  i/0  =  1.  Further,  o„,„  = 

±>/» -0-2-1-  The  sign  is  defined  by  the  condi¬ 
tion  detO  —  1.  The  quantities  i /,,  tx,-  are  called 
Schur  parameters  and  auxiliary  Schur  param¬ 
eters,  respectively.  As  we  saw  above  the  Schur 
parameters  i/j,  *  =  1,  1,  determine  0 

uniquely. 

On  the  other  hand,  if  we  know  the  entries 
=  ^i/^t+ii  *  =  1  •  **,n  —  1  of  the  matrix 
O  we  can  determine  7,  by  (2.5)  up  to  a  sign.  In 
other  words,  the  entries  o.+i,,-  (auxiliary  Schur 
parameters)  determine  O  almost  uniquely. 

3  Explicit  formulas  for  the 
evolution  of  auxiliary  Schur 
parameters  under  the  Toda 
flow 


Proof:  We 

know  [7]  that  0(t )  =  £(t)C?(0)fl(t)_1,  where 
exp(0(O)f)  =  Q(t)R(t),  Q(t)  is  orthogonal, 
and  R(t)  is  an  upper  triangular  matrix  with 
positive  entries  on  the  main  diagonal.  We  then 
clearly  have 


_  _  r«+l,*+l(0  _  m\ 


i  =  1, •••,«-  1.  Here  R(t)  =||  ry(<)  ||.  The 
operator  A*  R(t)  naturally  acts  on  the  i-  th 
exterior  power  f\'  Rn  by  the  following  rule: 
A'  R(t)(vi  A  ...  A  v.)  =  R(t)v  1  A  . . .  R(t)vi  for 
any  t?i,  •  •  • ,  €  Rn.  We  have,  further,  the  fol¬ 

lowing  relations: 


rjj(t)  =<  R(t)ei ,  R(t)ei  >= 


<  exp(O{0)t)ei,exp(O(Q)t)ei  >= 


Let  O(t)  =||  Oij(t)  ||  be  the  solution  to  the 
Toda  flow 

o  =  [o,xO], 

such  that  0(0)  is  upper  Hessenberg  orthogonal 
and  irreducible.  Here  irO  =  0_  -  0T_  and  0_ 
is  strictly  lower  triangular  part  of  O.  Then 
O(t)  possesses  the  same  properties  and  O(t) 
converges  when  t  — ►  00  to  a  block  diagonal 
matrix.  Each  two  by  two  block  corresponds  to 
a  pair  of  complex  conjugate  eigenvalues.  The 
blocks  are  arranged  in  the  decreasing  order  of 
real  parts  of  eigenvalues  [7,  5).  From  the  pre¬ 
vious  discussion  we  know  that  0(t)  is  almost 


<  ej,r(/)ei  >=  Ai(/). 

And  more  generally 

r?i(0  •••*$■(<)  = 


<  ex  A...Ae„/\r(*)(ex  A  ...  A  e,)  >= 
*  =  1,  •  •  • ,  n.  By  (3.2)  we  easily  obtain 


A.(0, 

(3.2) 


rj+U+i(t) 

*m(0 


The  result  now  follows  by  (3.1).  ■ 
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We  have  the  following  differential  equations  for 

Oi  : 

G\  —  l.t+l  O, 

*  =  1,  •  •  • ,  n  —  1.  Recalling  that  Oj,,  =  — 
i  =  1, •••,«,  i/n  =  ±1,  and  v?  +  a}  =  1,  we 
obtain  i/iU,  +  d,tr ^  =  0  or 

di  =  -  J'.+i), 

~  Vi+ 1), 

t  =  1,  •  ••  ,n  —  1.  It  is  interesting  to  find  how 
moments  rv(t)  =<  ei,0(t)‘ei  >  evolve  under 
the  Toda  flow.  Consider  the  family  of  rational 
functions 

ft(z)  =<  eu[zE  -  0(t))~lei  >  . 
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